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Abstract: We investigate the rates for multielectron recombination within a dense plasma environ- 
ment in local thermodynamic equilibrium (LTE). We find that these multielectron recombination rates 
can be high within dense plasmas, and they should be treated in the simulations of the plasmas created 
by intense radiation, in particular for plasmas created by intense VUV radiation from free-electron- 
laser (PEL) or for modelling the inertial confinement fusion (ICF) plasmas. 



1 Introduction 

Electronic many-body recombination processes that occur within dense, cold plasmas have attracted 
much attention recently as a possible explanation of the strong energy absorption within laser-created 
plasmas [1,2]. During the first cluster experiment performed at the free-electron-laser FLASH facility 
at DESY [3, 4] such plasmas were created after the irradiation of the xenon clusters with VUV photons 
of energy, = 12.7 eV. Pulse duration did not exceed 50 fs, and the maximal pulse intensity 
was, I < 10^^ W/cm^. Highly charged Xe ions (up to +8) of high kinetic energies were detected, 
indicating a strong energy absorption that could not be explained using the standard approaches [4- 
6]. More specifically, the energy absorbed was almost an order of magnitude larger than that one 
predicted with classical absorption models, and the ion charge states were much higher than those 
observed during the irradiation of isolated atoms at the similar conditions. This indicated that at 
these radiation wavelengths some processes specific to many-body systems were responsible for the 
enhanced energy absorption. 

The lowest order many-body recombination is three-body recombination, where two continuum 
electrons are involved initially, and in the final state one of these electrons is captured by the target ion. 
The excess energy released by the recombining electron is then carried away by the other outgoing 
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electron so that the three-body recombination does not involve any emission of photons [7]. For this 
to happen it is necessary that the electron density must be high enough so that the collisional processes 
are more probable than the dissipative radiation processes [8]. The three-body recombination is the 
inverse of the electron collisional ionization, and this relation is used in deriving the relevant cross- 
section formula, by detailed balance. 

In the similar way higher-order recombination processes can be described. More than two elec- 
trons are then involved. Rates for these processes strongly depend on the electron density, and can 
again be calculated with the rates of the corresponding inverse many-body collisional processes. 

The mechanism of plasma heating due to many-body recombination that was proposed in Ref. [1] 
in order to explain the experimental results of Ref. [3] can be best described as a sequence of Auger 
processes followed by photoionization. According to [1], the strongly coupled nanoplasma is created 
shortly after the exposure. Within this plasma, a probability of finding two or more electrons being 
close to an ion is relatively high. One electron can be captured by the ion, while the other electrons 
(spectators) carry away the excess energy released by the recombining electron. The captured electron 
may then absorb an radiation photon and become ionized once more. This cycle repeats many times 
during the duration of the pulse, leading to an efficient heating of the electron cloud. Ions of charge 
up to +7 were predicted with this model for the Xeiooo cluster. 

Here we aim to study a different many body recombination process that may also occur within 
dense plasmas. This is multielectron recombination (MER) that is an inverse process to multiple col- 
lisional ionization. Several electrons then recombine simultaneously, and a single spectator electron 
carries away the excess energy. We will evaluate the significance of this process within plasmas cre- 
ated by VUV FEL radiation from xenon and argon clusters [3] and also within inertial-confinement- 
fusion plasmas [9]. We will first derive analytical formulas for MER cross-section and rates, using 
the microscopic reciprocity relation and the detailed balance principle [10]. Using these formulas, we 
will obtain numerically the MER rates for Xe and Ar plasmas. These rates will give an estimate, how 
significantly these processes contribute to the dynamics within dense plasmas. 



2 Calculation of rates for multielectron recombination 

MER is an inverse of the multiple collisional process (Fig. I), where the primary electron ionizes an 
atom, releasing one or more secondary electrons. At n = 1 we arrive at the three-body recombination 
(one-electron recombination), at = /c we have k-electron recombination. In what follows we will 
use the notation as in Ref. [10]. 

Derivation of the MER cross sections and rates bases on the validity of the quantum mechanical 
reciprocity relation: 

w{i^ f)=w{f ^i), (1) 
where w{i — > /) is the transition probability per unit time for the transition \i >^ |/ > between 
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Figure 1 : Collisional processes within a dense plasma: (a) multiple collisional ionization and (b) its 
inverse process: multielectron recombination. The electron, e, is a spectator. 



single quantum states of the states i and /. Reciprocity relation (1) is an essential condition for 
the validity of the principle of detailed balance [10]. The conditions for validity of the reciprocity 
relation can be obtained within the framework of the general scattering matrix theory [10]. The result 
is that the reciprocity relation w{i /) is not generally vaUd, not even if the Hamiltionian of the 
system is invariant under space reflection (P) and time reversal (T). However, relation (1) holds in 
the framework of perturbation theory, and thus applies to radiative processes and to colMsions in the 
Born approximation. The reciprocity relation holds moreover for all collision processes involving 
unpolarized particles, provided that the Hamiltonian is P- and T-invariant. In what follows we will 
consider only unpolarized particles, therefore the reciprocity relation will be valid for the processes 
considered. 

Explicit form of the reciprocity relation for specific processes can be derived directly from quan- 
tum mechanics. Here, following Ref. [10], we will derive it through the principle of detailed balance, 
by writing down the rate equations and using the explicit form of the thermal distribution functions. 
This procedure is simple and straightforward, while the relations so found are independent of this 
method, and are identical with those derived directly from quantum mechanics. We will assume the 
electron distributions to be isotropic. 

Reaction equation for a collisional process (/c-body multiple ionization) and its inverse process 
(fe-body multielectron recombination) in Fig. 1 can be written as: 



^(^o) + e{E) ^ A+{E+) + e{E') + e(^i) + . . . + e{Ek 



(2) 



where A denotes atom or ion in the state of initial energy Eq, A+ denotes the ionized atom or ion in 
the state of energy £'+, e is electron of energy E{Ei). Energy of the spectator electron changes from 
E ^ E' during the process. This reaction leads to the balanced rate equation: 



Po\l — Pefe{E)no+{E;E',Ei, 

rUe 



■ dEk = 



P+ 



— pe fe{E') dE'J^ MEi) dEl 
me V rUe 



2E„ 



Pe U{Ek) dEk^+o{E', Ei,...,Ek;E) dE, 



(3) 



where we introduced the following notation: po denotes density of atoms/ions in initial energy state 
Eq, denotes density of atoms/ions in the final energy state pe is the electron density and rUe is 
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the electron mass. The distribution, fe{E) is a normaUzed isotropic Maxwell-Bohzmann distribution 
of electron kinetic energy, and quantities Qq^{E; E', Ei, . . . , En) and i7+o(-£'') Ei, . . . , En', E) are 
differential cross sections (in respect to energy) for the multiple coUisional ionization and multielec- 
tron recombination respectively. The angular distribution of electrons was assumed to be isotropic, 
and the integration over scattering angles was performed. 

The rate for multielectron recombination is defined as: 

dE' dEi ... dEk dE, (4) 
and can be related to the collisional multiionization rate: 

f I2E 

Rk(i) = Pe / \ — fe{E) no+{E; E',Ei,...,Ek) dE' dEi . . . dEk dE, (5) 
through a relation that follows from Eq. (3): 

P+ Rk{RC) = Po Rk{i) ■ (6) 

In a balanced state the relation between the atom/ion distributions in states £'0 and E^ is given by 
a Saha relation [10] known also as principle of detailed balance. For k-multielectron recombination 
and k-multiple ionization this relation can be derived (not shown here) from Eq. (2) as: 

k 



P+ 9+ V2 



where jS = 1/kT is the inverse of the temperature within the system and Ae = , ^ is the 
thermal de Broglie wavelength of free electrons. Coefficients, go,g+, are statistical weights of the 
bound energy levels, Eq,E^, respectively, and £'+0 = ^+ ~ Eq is the energy needed for transition 
from state Eq to state £+ (ionization energy). After substituting Eq. (7) to Eq. (6), we obtain: 



RkiRC)=Rkii)f^ef'^+° (^^PeXl) . 



k 

(8) 



9+ 

From Eq. (3) also a microscopic relation between ionization and recombination cross sections can 
be obtained: 

24fc 

goEno+{E;E',Ei,...,Ek)=g+j^{7rme)''E'Ei . . . E^ n+o{E' , Ei, . . . , Ek; E). (9) 

This is a generalization of the Fowler relation obtained for three-body recombination [10, 11]. As 
expected, at A; = 1 Eq. (9) reduces to the Fowler relation: 

gQEno+iE;E',Ei)=g+-^{Trme)E'Ein+o{E',Er,E). (10) 

For plasmas in a non-equilibrium state (non-LTE plasmas) Eq. (8) is not longer valid. The non- 
equilibrium rate for multielectron recombination can then be still obtained from Eq. (4), using the 
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recombination cross section derived from the microscopic relation, Eq. (9), and convoluted with 
non-equilibrium electron distributions. The difficulty lies in the correct parametrization of the differ- 
ential multiple ionization cross section, which is not known for an arbitrary value of n > 1 [12, 13]. 
Therefore we restrict here to LTE plasmas, for which the direct rate equation, Eq. (8), can be used. 

3 Estimation of multielectron recombination rates for xenon and argon 

We used Eq. (8) to estimate the rate for multielectron recombination within plasmas in thermal equi- 
librium. In order to calculate the collisional ionization rate, we used experimentally determined total 
cross sections for collisional multiionization [14]. 

Eq. (8) implies that the electrons recombine in most cases to the higher Rydberg states, for which 
the statistical weights, go, are large. Therefore here we estimate the recombination rates to an "aver- 
age" Rydberg state. This is the state closest to the average of the available states within the plasma. 
This average is obtained with weights that correspond to the statistical weights, go, of the considered 
states. The recombination rates obtained are plotted in Fig. 2. Density range, p = 10^^ — 10^"^ cm~^ 
at temperatures T = 1 — 10 eV corresponds to plasmas created by intense VUV laser pulse from EEL 
[15-18]. The rates obtained for densities p = 10^^ — 10^^ cm~^ and temperatures T ~ 100 eV are 
relevant for inertial confinement fusion plasmas. 

At the density, p = 10^^ cm~^, and temperatures, T = 1 — 10 eV, the highest rates are: the 
rate for 3-electron recombination (Xe), and 2-electron recombination (Ar). This tendency keeps with 
increasing electron density. The absolute values of recombination rates (ra = 1 — 3) at p = 10^^ cm~^ 
are then between 10~^ — 10^ 1/fs for Xe and 10~^ — 10^ 1/fs for Ar, depending on the temperature. 

In the region of high densities, p > 10^"^ cm~^, and at high temperatures, T = 100 eV, relevant 
for ICE plasmas, the absolute recombination rates are high: 10^ — 10"^ 1/fs for Xe, and 10^ — 10^ 1/fs 
for Ar at p = lO^^ cm The I -electron recombination is predominant but the other recombination 
rates (n = 2, 3) are less by the factor of 10 (Xe) and 10-100 (Ar). 

4 Summary 

The aim of our study was to estimate the multielectron recombination rates (ra > 1) within dense 
plasmas. 

We derived the MER recombination rate for LTE plasmas built of xenon or argon atoms, using 
the experimental cross sections for multiple ionization of atoms. We considered the recombination to 
the average Rydberg state. Plasma effects as lowering of ionization thresholds, modifications of the 
cross sections have not been treated. Cut-off of available atomic states due to the plasma environment 
was included. 



5 




o 
ir 



1e-20 
1e+20 



1e+21 1e+22 

Pe [cm 



1e+23 

3i 



1e+24 1e+25 



(b) 




o 
ir 
CC 



1e-20 
1e+20 



(c) 



1e+21 1e+22 1e+23 1e+24 1e+25 
Pe [cm"^] 



1e+08 
1e+06 
10000 
100 
1 

0.01 
1e-04 
1e-06 
1e-08 
1e-10 
1e+22 



, T=1 eV 


Xenon 


1 T=10 eV 




■ T=100 eV 


















T-el. 




2- el. 

3- el. ■ ' 





1e+08 
1e+06 
10000 
100 

1 

0.01 
1e-04 
1e-06 
1e-08 
1e-10 
1e-12 
1e+22 



1e+23 
Pe [cm"^] 



(d) 



1e+23 
Pe [cm"^] 



1e+24 



, T=1 eV 


Argon 


1 T=10 eV 




T T=100 eV 












— ' 7.— s'""*^"' 










^ T-el. — 




2- el. 

3- el. ■ - 





1e+24 



Figure 2: Rates for multielectron recombination: 1, 2, 3-electron recombination to average Rydberg 
state as functions of plasma density for plasma in LTE. They are plotted for: (a) xenon with a closeup 
(b), and (c) argon with a closeup (d). Arrows show how the rates change with increasing values of 
electron temperature, T = 1, 10, 100 eV. 
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We showed that in the p ~ T region relevant for plasmas created after irradiation of Xe and 
Ar clusters with intense VUV pulses from FEL, MER recombination dominates over the 1 -electron 
recombination. The values of recombination rates (n = 1 — 3) are between 10^^ — 10^ 1/fs for Xe 
and 10"*^ — 10^ 1/fs for Ar at p = 10^^ cm~^, depending on the temperature. In the p — T region of 
ICF plasmas, the absolute values of MER are high (10^ - 10^ 1/fs for Xe, and 10^ - 10^ 1/fs for Ar 
at p = 10^^ cm~^). The l-electron recombination dominates in this regime but the 2- and 3-electron 
recombination rates are less by factor of 10-100. 

Our estimations imply that the MER rates should be treated in simulations of processes occurring 
within ICF plasmas and FEL laser created plasmas. The quantitative validation how significantly the 
MER processes contribute to the plasma dynamics would require a detailed rate equation approach, 
where both competing processes: non-sequential ionization and MER recombination are included, 
and plasma effects have to be treated. This approach is, however, beyond the scope of the present 
study. 
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